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FOREWORD 


Today petrologists, mineralogists, and cry stallograp hers are the possessors of very-well-equipped 
laboratories. The laser probe, the electron microprobe, and the ion microprobe mass analyzer can pro- 
duce results within minutes, results that used to take days and weeks to attain. Besides X-ray methods, 
there are several spectroscopic techniques, such as Mossbauer spectroscopy, that can be used to exam- 
ine the distribution of cations over nonequivalent structural sites in a crystal. The experimental 
methods of synthesizing mineral assemblages in the laboratories have been greatly improved. Refined 
techniques exist for the control of gas fugacities and the creation of low to very high pressure and 
temperature conditions in petrological experiments. These technical capabilities are further enhanced 
by the use of computers, which can analyze numerical data and the consistency or inconsistency of 
experiments. 

As a result of increased experimental capability, phase equilibria data have been gathered both 
from efforts to synthesize mineral assemblages in the laboratory and from chemical analysis of natural 
mineral assemblages. To make a meaningful analysis of both these types of data, thermodynamic data 
on rock-forming phases and crystalline solutions must be available. Unfortunately, obtaining these 
quantities by thermo chemical and calorimetric methods has been a very slow process, and the possi- 
bility of obtaining such quantities by other methods must be considered. 

In suitable cases, retrieving thermodynamic data from the phase diagrams themselves may be 
considered. These quantities may be considered significant if they are obtained consistently from dif- 
ferent phase diagrams. However, in the case of crystalline solutions, thermodynamic functions of 
mixing cannot be obtained without the use of certain solution models. Two types of phase diagrams 
may be considered. The first type is the diagram with the solvus or binodal curve bounding a -binary 
two-phase region. The analytical methods have been discussed by Thompson (1967), Thompson and 
Waldbaum (1969 a,b), and Green (1970). These methods, however, are of limited use for rock-forming 
silicates because many crystalline solutions do not show any solvus relationship. The second type of 
phase diagram is the Roozeboom-type figure in which the concentration of a component in one phase 
is plotted against the concentration of the same component in the coexisting phase. Such a distribu- 
tion relationship, based on simple ion-exchange reactions, was discussed initially by Ramberg and 
DeVore (1951) and later by Kretz (1959) and Mueller (1960). It is possible to retrieve useful thermo- 
dynamic information from the distribution data in ion exchange collected by Nafziger and Muan 
(1967), Larimer (1968), Medaris (1969), and Schulien et al. (1970), among others. 

Useful thermodynamic information may also be obtained by considering distribution of cations 
within the crystalline solution. Ghose (1961) found an interesting Fe 2+ -Mg 2+ distribution in 
cummingtonite. Since then, Fe 2+ -Mg 2+ order-disorder has been studied in several silicates by crystal- 
lographers. Theoretical framework for considering the homogeneous equilibria of intracrystalline 
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cation distribution has been presented in papers by Mueller ( 1 962), Matsui and Banno (1965), 

Perchuk and Ryabchikov (1968), Thompson (1969), and Grover and Orville (1969). 

The purpose of this work is to discuss these methods of obtaining thermodynamic quantities 
and some aspects of partitioning elements in coexisting phases by considering the definition of crystal- 
line solutions, the definition of components in a silicate mineral, and the definition of chemical 
potentials of these components. The solution models involved are also considered. An example of cal- 
culating thermodynamic functions of mixing in the CaW0 4 -SrW0 4 system is given. 

It is hoped that this work will generate enough interest among fellow scientists to gather useful 
ion-exchange data on coexisting phases and among crystallographers to gather data on site occupancies 
in rock-forming silicates. 
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SYMBOLS 


Superscripts are generally abbreviated names of the minerals to which the thermodynamic func- 
tions are ascribed. Subscripts refer to components of the crystalline solution or the chemical system. 

a\ activity of component A in phase a 

A Q , A v A 2 energy constants in equation for excess free energy of mixing expressed as a polyno- 
mial in mole fraction 

A, B, C used as a subscript denotes components A, B, and C 

A, B, C energy constants used in equations describing the relation between activity and mole 

fractions 

/ activity coefficient 

G molar Gibbs free energy 

G £M excess molar Gibbs free energy of mixing 

G, m ideal molar Gibbs free energy of mixing (= Z) R T In xj 

G u total molar Gibbs free energy of mixing 

H molar enthalpy 

i/ EM excess molar enthalpy of mixing 

H u molar enthalpy of mixing 

K thermodynamic equilibrium constant 

K d distribution coefficient 

Ml, M2 structural sites in the crystal 

Ti Avogadro’s number 

P pressure 

q contact factor 

R gas constant 

S molar entropy 

^EM excess molar entropy of mixing 
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T 

T 

C 

W 

w 

W' 

X ? 

1 

z 

M AM 

aim 

bi 

en 

gar 

ol 

opx 


ideal molar entropy of mixing 
molar entropy of mixing 
absolute temperature 
critical temperature 

energy constant or the interchange energy used in the regular solution model 

ft w where w is a function of P and T as in the simple mixture model 

ft w where w is independent of P and T 

mole fraction of a component i in phase a. 

coordination number 

chemical potential of a component i in phase a 
chemical potential of a pure component AM 
almandine: Fe 3 Al 2 Si 3 0 12 
biotite: K(Fe, Mg) 3 AlSi 3 O 10 (OH) 2 
enstatite: MgSi0 3 or MgMgSi 2 0 6 
garnet: (Mg, Fe, Ca, Mii) 3 Al 2 Si 3 0 12 
olivine: (Fe, Mg) 2 Si0 4 

orthopyroxene: (Mg, Fe)Si0 3 or (Mg, Fe) 2 Si 2 0 6 
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Chapter 1 


THERMODYNAMIC RELATIONS IN CRYSTALLINE SOLUTIONS 


Thermodynamic relations between the concentration of a component in a solution and its chemi- 
cal potential and other thermodynamic functions of mixing are presented here. The details of the 
simplifying assumptions and the methods of statistical thermodynamics have been given by Denbigh 
(1966), Guggenheim (1952, 1967), and Prigogine and Defay (1954), among others. Recently Thomp- 
son (1967) also considered the properties of simple solutions. Besides a summary of thermodynamic 
relations in solutions, the difficulties encountered in their application to silicate minerals will be con- 
sidered. Some of these problems, such as the choice of a component and definition of its chemical 
potential in a silicate, have been discussed by Ramberg (1952a, 1963), Kretz (1961), and Thompson 
(1969). 

CRYSTALLINE SOLUTIONS 

The crystalline solutions considered here are rock-forming silicates forming isomorphous series 
with one another. Such crystalline solutions have a definite structural framework with generally two 
or more kinds of nonequivalent structural sites. The type of sites and the ions that occupy them vary 
in different crystalline solutions. The overall crystal symmetry of a solution does not change as a func- 
tion of the composition, though certain microscopic details within the crystal, i.e., the form and size 
of the individual structural sites, may change with changing composition. 

Orthopyroxene (opx) (Mg, Fe) 2 Si 2 O g may be considered as an example. In the crystal structure 
there are single silicate chains parallel to the c-axis held together by the octahedrally coordinated Mg 2-1 " 
and Fe 2+ . There are two kinds of structurally nonequivalent sites M 1 and M2 occupied by Mg 2+ and 
Fe 2+ . The Ml octahedral space is nearly regular polyhedral, but the M2 space is quite distorted. As a 
result of varying Mg 2+ and Fe 2+ in the composition of the crystal, the general symmetry of the crystal 
does not change, but there are distinct changes in M 1 and M2 polyhedra. The former becomes more 
regular and the latter more distorted with increasing Fe/Mg ratio. Such microscopic changes at the 
structural sites within the same crystal framework may be regarded as continuous; and the resulting 
energy changes, a consequence of the mixing or solution of the species to form a crystalline solution. 

CHOICE OF A CHEMICAL COMPONENT 

The definition of a component in a mineral is not unique. The components in orthopyroxene may 
be considered to be the molecules MgSi0 3 and FeSi0 3 or the molecules MgO, FeO, and Si0 2 or the 
ions Mg 2+ , Fe 2+ , Si 4+ , and O 2- . In petrological studies, the choice of a component is determined by 
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known or postulated chemical reactions involving a mineral. In such studies, the use of components 
such as FeSi0 3 or FeO is convenient, even though there are no discrete units of this kind in the ortho- 
pyroxene crystal structure. However, when the thermodynamic properties of silicate crystalline solu- 
tions are being considered, it is only realistic to consider the ions as the components. (See Bradley, 
1962.) Indeed, it can be noted that if the substitution of the cation Fe 2+ by Mg 2+ in orthopyroxene 
does not produce any changes in the silicate framework or if there are any slight changes, they are 
directly a function of the changing Fe/Mg ratio; the alternative methods of defining FeSi0 3 or Fe 2+ as 
a component are equivalent. (See also Saxena and Ghose, 1971.) 

CHEMICAL POTENTIAL AND ACTIVITY OF A COMPONENT IN A MINERAL 

A solution is ideal if the chemical potential of every component is a linear function of the loga- 
rithm of its mole fraction according to the relation 

/J- i =/U- 1 + R T lnx ; (1-1) 

where fx i is the chemical potential of i in a solution and p* is the chemical potential of pure i. p 1 is a 
function of pressure (P) and temperature (T) only. In a binary solution a whose composition is 
(A, B)M, where M represents the anion group or the silicate framework and A and B, the cations that 
substitute for each other, there is a choice between adopting the cations A and B as components or the 
end member molecules AM and BM. As noted before, under certain conditions, the mole fractions 
may be calculated as 

x a = A 

A aTb 

or 

„a AM 

X m AM + B M ' 

These expressions could be considered equivalent to each other. For chemical potentials, 

Pa + R T ln *A 0-2) 

or 

RTl n*M’ 0 - 3 > 

where p AM and ai am are chemical potentials of A and AM in a standard state. The standard state AM 
is well defined, but the standard state with reference to cation A needs definition. In orthopyroxene, 
this is like referring to the chemical potential of Mg 2+ in pure (Mg, Mg)Si 2 O g . The Gibbs free energy 
for the pure end member MgSi0 3 is defined and experimentally measurable, but the meaning of free 
energy of Mg 2+ in pure enstatite (en) is little understood and experimental methods remain to be 
developed for its measurement. 

In theoretical discussions, however, in which the measured values of the potentials are not of con- 
cern, the definition of chemical potential of a cation in a crystalline solution is not only permissible 
but also useful. Kretz (1961) defines the chemical potential of Mg in orthopyroxene as 
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Ms Ip, T, n , n , n 
Fe Si 0 


(1-4) 


where n is the number of cations in the formula. 

In many crystalline solutions, when their compositions are expressed in the simplest form, there 
are two or more cations in one mole. Examples are olivine (ol) [(Fe, Mg) 2 Si0 4 ] and garnet (gar) 
[(Fe, Mg) 3 Al 2 Si 3 0 12 ] . The chemical potential of a component using the molecular model is ex- 
pressed as 


and 


ol Fe 2 Si0 4 ol 

^Fe 2 Si0 4 = ^ + R T ln *Fe 2 Si0 4 


(1-5) 


ngar _ a Fe 3 A1 2 Sl 3°i 2 n T , gar 

Pe 3 Al 2 Si 3°12 ^ RT F e 3 A1 2 Si 3 0 12’ 


( 1 - 6 ) 


where the first p in both equations is the chemical potential of the end member in the solution and the 
second p is the chemical potential of the pure end member. If the cation Fe 2+ is considered as a 
component, 


M 


ol 

Fe 


= M 


Fe 2 Si0 4 

Fe 


+ 2 R T In X° l 


(1-7) 


and 


„„„ Fe,Al -Si ,0, , 

Mg“ r =M F e 3 2 3 12 + 3 RT ln*|r , 


( 1 - 8 ) 


where the first p in both equations is the chemical potential of Fe 2+ in the crystalline solution and the 
second p is the chemical potential of Fe 2+ in the pure end member. The mole fractions x are the same 
quantities in both the molecular and ionic models. It may be desirable to consider the chemical for- 
mula on a one-cation basis; i.e., olivine is considered to be (Fe, Mg)Si 0 5 0 2 and garnet to be 
(Mg, Fe)Al 2 / 3 Si0 4 . In this case, 


or 


and 


or 


M F ; = 2M 


o 1 


MFe Si, 5 o 2 = 2 M 


M 


Fe ^ 


/X% ar = Lll 

FeAl 2 / 3 Si0 4 3 M 


Fe 2 Si ° 4 +/?rln ,o. 

(1-9) 

: Fe2Si ° 4 + R T lnx°l Si 0 

0.5 2 

(1-10) 

,Al 2 Si 3 0l 2 + R T In*!*' 

r c 

(1-11) 

, Al ,Si ,0 , , 

3 2 3 1 2 + R T In *ear 

FeAl SiO 

(1-12) 
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The usefulness of these relations is mentioned later in connection with the composition of 
coexisting minerals. 
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The activities of the components will be the primary concern of this document. For a binary 
ideal solution, the activity is equal to its mole fraction. In olivine the activity of the fayalite (fa) mole- 
cule is 

«fa = 0-13) 

or for Fe 2+ , 

. ,.v : ; ' % e = (<) 2 - (1-14) 

Similarly for garnet and almandine (aim), 

a i 
aim 

and 

« Fe = 0|D 3 • 0-16) 

It is desirable to consider many reactions, particularly the ion-exchange reaction, on a one-cation 
basis; i.e., to consider olivine as (Fe, Mg)Si Q 5 0 2 , etc. Activity of a cation is then equal to its mole 
fraction. It is necessary to specify that although in this situation x Fe [the mole fraction Fe 2+ /(Fe 2+ + 
Mg 2+ )] is numerically the same as x fa (the percent of fayalite), the activities are different. x Fe = (a Fe ) 1/2 
in the ionic model, but jc fa = a fa (the activity of fayalite) in the solution. 

NONIDEAL BINARY SOLUTIONS 

The relation between the chemical potential of a component i and its activity in a solution is given 
by 

t /j . j .= /u 1 + R T In a.. (1-17) 

The ideal solution is the limiting case when a i is equal to the mole fraction x { . In all other cases, the 
relation between a y and x i may be expressed as 

a . = f . x . , (1-18) 

where /. is the activity coefficient of the component i in the solution. 

The free energy of mixing G M for a binary solution (A, B)M is given by 

C m = *a R T ln a A + *b R T ln °B 

= R T (* A ln X \ ln + R T (*A ln f A + X B ln f B> 

= C im + 6 em . (1-19) 

The first term, G IM , is the ideal free energy of mixing, and the second term, G £M , is the excess 
free energy of mixing that is due to the nonideal nature of the system. G EM is one of the functions of 
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mixing termed ‘excess functions.’ Details of the excess functions of mixing may be found in Prigogine 
and Defay (1954) and Thompson (1967). 

REGULAR SOLUTION MODEL 

The thermodynamic properties of crystalline solutions of mineral assemblages either obtained in 
experiments or occurring in rocks must often be calculated. In many situations the use of certain 
models for the activity composition relationship helps to assess such properties closely. The regular 
solution model of Guggenheim (1952) is next in simplicity to ideal solution model. 


Zeroth Approximation 


The excess free energy of mixing G £M in a regular solution with the zeroth approximation, i.e., 
the approximation of complete disorder, is given by 

G EM=*A*B r ’ d' 20 .) 

where A and B are components of a solution (A, B)M, W' is equal to Tl w, and Tl is Avogadro’s number. 
W' is often referred to as the interchange energy. Regular solutions are very important in this work; 
therefore, the parameter w will be briefly discussed. A simplified account of this parameter is pre- 
sented by Denbigh (1966). It is assumed that the cations A and B are of roughly the same size and can 
be interchanged between lattice sites without change of lattice structure and without change in the 
lattice vibrations. There is interaction between A and B, given by the energy w, which is given by 


w 


= 2 uj ■ - • - 

AB AA 


BB 


d-21) 


where w AA is the increase in potential energy when a pair of A ions are brought together from infinite 
distance to their equilibrium separation in the solution. w AB and w BB are similarly defined. In spite 
of the interaction energy, it is assumed that the mixing of A and B is random. This means that the 
entropy of mixing is the same as that for an ideal solution and deviations are expressed entirely in 
terms of the heat of mixing. 

The thermodynamic equations for the regular solution model with zeroth approximation are 


and 


G EM ” ^EM 


= X.xJ," 

A B 


( 1 - 22 ) 



(1-23) 


The interchange energy W' is independent of P and T. Because the excess entropy of mixing is zero 
according to this model, the predictions of the values of and the heat of mixing H EM , which may 
often be different from G EM , are not satisfactory. 

Simple Mixture Model 

In the regular solution model W' is supposed to be independent of T and P. In Guggenheim’s 
(1967) latest version of the lattice theory, W' may be treated as an adjustable constant required to fit 
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the experimental data to the model. Such an energy parameter with a symbol W may be called a 
cooperative free energy. 2W is in a sense the free energy increase in the whole system when an A-A 
pair and a B-B pair are converted into two A-B pairs. It is expected that if W is fitted to the free 
energy data at each temperature, the large errors usually found in the predictions of <? EM and H EM 
with composition may be at least partly eliminated. For a random mixing approximation, the various 
excess functions are given by 


G EM ~ *A*B ^ ’ 


c 3 W 

— S w = x A X n 

EM A B j j 1 


and 


H EM = *A*B ~ T 




The activity coefficient is related to the mole fraction by 

, , If 2 

In / . = - — 

A R T B 


(1-24) 

(1-25) 


(1-26) 


(1-27) 


Quasi-Chemical Model 


The main assumptions required for this model are similar to those of the regular solution model 
in the preceding sections. Only the configurational partition function of the solution contributes to 
the thermodynamic excess functions. The intermolecular forces are central and short range, and there- 
fore the internal energy at 0 K may be obtained by an addition of the pair potentials. The assumption 
of complete randomness is not required here. Therefore any differences found in the calculated values 
of the excess functions of mixing by the zeroth approximation and by the quasi-chemical approxima- 
tion are the result of ordering considered in the latter. 


In binary solutions for which the two components A and B are of similar size, the activity coeffi- 
cients are given by the equations 




/3 + 1 - 2 x ] 

x A O 3 + 1) 




(1-28) 


and 


z 


f B 


P ~ 1 + 2 X B 

* A 0 + 1 ) 


(1-29) 


where z is the coordination number and (3 is given by 

/3 = [i + 4x A x B (e 2W/zRT - 1)] 1/2 . 


(1-30) 
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/3 = 1 for a perfectly random mixture. If (3 > 1, a tendency for clustering exists, and if /3 < 1, a trend for 
compound formation exists. 


G em andi/ EM are given by 


C - — z R T 
EM .2 


/3+1-2*b fi 1 + 

x A In - - — + * B In- 


A * A 0 + 1) 


*A 0 + J ) 


and 


^ EM ~ p + i-* A * B 


('-If)- 


(1-31) 


(1-32) 


The various equations of the quasi-chemical approximation may be expanded as power series in 
2W/zRT: 


G m = 


"l- 1 

( ^ W ) 

1 / 2 W V 

-v r — I 1 

2 

V R t) 

A b 6 \z /? ry 


K ' v r) 2 + 


(1-33) 


and 


f* = w *l 


i + I 

/2 r\ 

2 

V** t ) 


l 2 W 


v V 


6 V z R T 


"a (1-30 hr; *a (l-H* B + 28*a-20*3 ) + ... 


(1-34) 


/ B may be obtained by replacing A by B in Equation 1-34. 


For molecules that are not very similar in size, a contact factor must be included (Guggenheim, 
1952, p. 216) in these equations to take the size differences into account. The contact factors may be 
found roughly proportional to the molar volumes or ionic radii. The activity coefficients are given by 


u = 


1 + 


<£ b 0 - !) 

4>. 0 + 1) 


*9 «/2 


(1-35) 


and 


f B = 


1 + 


0 A 0- !) 

4> n 0 + 1) 


z Q „/2 


(1-36) 


where q A and q R are contact factors related to the contact fractions <p A and 0 B and the mole fractions 
x A andx B by 


= 


*A?A 

X A 9A + X B ?B 


(1-3 7a) 
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and 




B “B 


C A 9 A + *B?B 


( 1 -3 7 b) 


For more details on the derivation and significance of the constants <? A and q h and the fractions <j> A 
and 0 B , reference may be made to Guggenheim (1952, p. 216) and King (1969, p. 488). 0 in Equations 
1-35 and 1-36 is obtained by replacing x A and * B by <j> A and 0 B , respectively, in Equation 1-30. 


The other excess functions are given by 


C EM = £ » *T 


c a9a 


In 


1 + 


^ (P ~ 1) 

4> k 08 + 1 ) 


+ *6 98 ln 


1 + 


(P ~ 1) 

(P + 1) 


(1-38) 


and 


H = 2 % _± * B ?A ?B ( r _ T B ^ 

(*A U + *B 9 B> + J ) \ 


3 T 


(1-39) 


5 em can be obtained by the standard equation 

^EM = ^EM “ ^ ^EM' 

SOME OTHER METHODS FOR REPRESENTING THE ACTIVITY-COMPOSITION 
RELATIONS 


In the case of a binary system of components A and B, the Gibbs-Duhem equation is 

x A d In f A + x B d In f B = 0 • (1-40) 

The changes d In f A and d In / B when due to composition change dx A at constant temperature may be 
written 


3 In /. 


Inf, 


+ x 


B\ 3 


0 . 


(1-41) 


A solution to the above equation was proposed by Margules in the form of a power series: 


and 


Inf. = a 


A*B + b A X l + C A *B + J A*B + 


(1-42) 


ln f B = fl B *A + , b B X l + C B X l + X X + • • • 


(1-43) - 
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When the series is terminated at x 3 , the following relations exist between the coefficients: 


and 


°B = 0 ’ 


b B “ b A + 3 C A / ' 2 + 2 d A + • ‘ ‘ ’ 


C_ (c, + 8 d:/3 + • • •) 


(1-44) 


Using these relations, Carlson and Colburn (1942) expressed the activity coefficients by the equations 

log/ A = (2 B - A) (1 - * A ) 2 + 2 (A - B) (1 -* a )3 (1-45) 


and 


log / B = (2 A- B) x\ + 2(B - A) . 


(1-46) 


Relations similar to these have been used by Thompson (1967) and Thompson and Waldbaum (1968, 
1 969a, b). 


Another two-constant equation is due to van Laar. The equation resulted from a theory based on 
the van der Waals equation of state. This theory is probably incorrect, but van Laar’s equation con- 
tinues to be useful for representing the activity-composition relation. This equation is 



(1-47) 


Similarly for the other component, 


lo *f B 



(1-48) 


For many chemical systems, van Laar’s equation provides a better representation of the data than the 
Margules two-constant equation. The relative merits of these two equations were discussed by Carlson 
and Colburn (1942). Finally, it may be remarked that a power series expansion as (see Equation 1-49) 
for the G em is now widely preferred. Therefore, only such expressions and not the equations men- 
tioned in this section will be used. Expressing C £M as a power series is a means of giving empirical 
description to deviations from the ideal, which is a better alternative to the power series expansions 
referring to individual activity coefficients mentioned above. (? EM expressed as a power series can be 
related more conveniently to other global properties of the mixture, such as the heat and volume 
change of mixing, than can the individual activity coefficients, which represent the deviations divided 
up, as it were, among the components. 
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GENERAL RELATIONS FOR BINARY CRYSTALLINE SOLUTIONS 


Excess functions in nonideal solutions may conveniently be expressed by a power series in the 
mole fraction. Guggenheim (1937) suggested that G EM can be expressed as a polynomial in x as 

C EM = * A l A 0 + A 1 ( X A ~ * B ) + Aj 1 ( *A ~ X B^ + * ’ *] * (1 ' 49) 

where A 1 , and A 2 are constants. When odd terms in Equation 1-49 vanish, the solution becomes 
symmetric. If A 2 and other higher terms are also zero, the simple mixture model with A () as the 
energy constant W in Equation 1-24 results. The expressions for the activity coefficients are obtained 
from 

<5 G fm 

RT ^f a = C em+ x bT ^L 

A 

= x* [A Q + a x (3 x k - * B ) + A 2 (x a - * B ) (5 x A * B ) + •.*•] . d-50) 

and 

3- C™, 

RTlnf B = G m -* A —- . 

° X B 

= X A K “ A 1 ( 3 X B ~ X A> + A 2 ( X B ~ X A> ( 5 * B _ * A) + ' • ’ (1 " 51) 

Equations for other excess functions of mixing may be derived from Equation 1-49: 


~ 5 EM - X A X B 


A Q I'd A 2 

3 ^nW‘*-‘»> + V 3 F 


( X A ~ x n) 2 + * 


(1-52) 


and 


^em ~ X A X B 


A o - T 


f * A C 

3 7 ' 


A x -T\ 


f d A i 


( X A - * B ) + 




f 3 A) 

+ 

A 2- 7 



( X A ~ X B ) 2 + • • 


(1-53) 
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Chapter 2 


THERMODYNAMIC STABILITY OF A SOLUTION 


INTRINSIC AND EXTRINSIC STABILITY 

A crystalline solution in an ideal state adds a certain amount of negative free energy of mixing to 
free energy of the system. With increasing positive deviations from the ideal state, this contribution 
becomes less and less. Below a certain critical temperature of unmixing, the solution unmixes to form 
two or more solutions. These energy changes obviously affect the stability of the entire system of the 
mineral assemblage. This instability of a crystalline solution, which is the result of the positive excess 
free energy of mixing, may be termed “intrinsic instability” (see Mueller, 1964). Ideal solutions are 
always intrinsically stable. A crystalline solution may also become unstable if the physical and chem- 
ical conditions change in such a way that certain reaction products form a lower free energy assemblage 
than the crystalline solution. This instability can be called extrinsic. A solution may be both intrin- 
sically and extrinsically unstable. The division is essentially artificial; however, it helps in understand- 
ing and describing certain petrological reactions as shown by Mueller ( 1 964). 

Olivine [(Fe, Mg) 2 Si0 4 ] and pyroxene [(Fe, Mg)Si0 3 ] may be considered as ideal binary solu- 
tions at high temperatures (~1400 K). In spite of their ideal character, orthopyroxenes with more 
than 5 5 mole percent of ferrosilite were found unstable at liquid-state temperatures by Bowen and 
Schairer (1935). The iron-rich pyroxene is unstable because of the instability of ferrosilite relative to 
fayalite and quartz. This is extrinsic instability. 

At low temperatures (~900 K) the situation is little different. Orthopyroxene is somewhat non- 
ideal, and high values of G EM are associated with the high ferrosilite content of the solution. The 
extrinsic instability of the solution relative to olivine and quartz is less because iron-rich pyroxenes 
(about 86 percent FeSi0 3 ) are stable in metamorphic rocks. The instability of pyroxenes with higher 
ferrosilite in metamorphic rocks may be due to both the extrinsic and intrinsic instability of the ortho- 
pyroxene solution. 

CRITICAL MIXING 

General Conditions 

The conditions for critical mixing in terms of free energy of mixing G M and the mole fraction x 
are 

d 2 C M /d x 2 - 0 ' (2-1) 
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and 

B 3 C u /d x 3 = 0. 

These may be expressed in terms of G EM as 


B 2 G m /^ x 2 = - R T/x (1 - x) 

and 

B 3 G^/d x 3 = - R T (2 x - 1) A 2 (1 - x) 2 . 


Simple Mixture 

For a simple mixture, 

where 


Gem = * (1 " *) V , 


= V (7*. 


P)- 


By successive differentiation of Equation 2-5, 

B 2 G^/B * 2 = - 2 w 

and 

B 3 C^/B x 3 = 0. 

By substituting Equations 2-3 and 2-4 into Equations 2-6 and 2-7, respectively, 

- 2 W = - R T/x (1 - *) 


and 

0 = R T (2 * - 1)A 2 (1 - *) 2 . 


These give the critical composition when x = 0.5 and 2 RT ( = W. 

General Nonideal Solution 

For a binary solution that is not a symmetric solution, 

C EM = * c 1 - *) po + A i (1 - 2 *) + A 2 (1 - 2 x) 2 + . . .] . 
Successive differentiation of Equation 2-10 with respect to x gives 

- — ™ = - 2 /1 0 - 6 A x (2 * - 1) - A 2 [10 - 48 * (1 - *)] 


( 2 - 2 ) 

(2-3) 

(2-4) 

(2-5) 

( 2 - 6 ) 

(2-7) 

( 2 - 8 ) 

(2-9) 

( 2 - 10 ) 

( 2 - 11 ) 
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and 


33 g em 

_= -12 A. + 48 4, (1 - 2x) . (2-12) 

Substitution of Equations 2-3 and 2-4 into Equations 2-1 1 and 2-12, respectively, gives equations that 
are transcendental and cannot be solved without a computer program using an iteration method. 


Formation of Miscibility Gaps in a Ternary Simple Mixture 

Consider a ternary simple mixture with components 1, 2, and 3. W for the three binary systems 
are W l2 , W y3 , and W 23 . The chemical potentials of the components in the solution are given by 


Mi = Mi (T, P) + R T In % x + R T In f 1 . . . , (2-13) 

where RT In / may be expanded in terms of x and W as follows: 

R T In f t = (x 2 ) 2 V 12 + (* 3 ) 2 V l3 + * 2 * 3 (IT 12 - W 23 + V u ) , (2-1 4a) 

RTinf 2 = (* 3 ) 2 y 23 + (* x ) 2 r 12 + * 3 (i^ 23 - r 13 + r la ) , (2-i4b) 

and 

« T In f 3 = (* t ) 2 * 13 + (x ,) 2 ir 23 + * 2 (Fj 3 - f 12 + F 23 ) . (2- 14c) 

At equilibrium in the two separated coexisting phases a and 0, 

M? (* 2 ’ * 3 . T) -m? (**, * 5 , T) = 0 . (2-15) 

H 2 and p 3 are similarly defined. Substituting Equations 2-13 and 2-14 into Equation 2-15 and rearrang- 
ing (see Kaufman and Bernstein, 1970, p. 226), 

R T In (*{/*?) + F la [(*$) 2 -(*“)] + F J 3 [(*f) -(*“) 2 ] + A I f (*5 *“) = 0, (2-1 6a) 

« r In (x^/x“) + F 12 [(1 - x£) 2 r (1 - *£) 2 J + ^ 13 [(*f) 2 -(*“) 2 ] 


- A W (1 - xP) - (1 - *“)] = 0 , (2- 16b) 


and 


R T In (*5/*S) + ^ 12 [(^) 2 - (*P 2 ]+ F 13 [(1 - *f) 2 - (1 - x“) 2 ] 

- A f (1 - x^) - x“ (1 - xpj = 0 . 


(2- 16c) 


where AW = W X2 + W 13 ~ ^23' 
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With the help of Equations 2-16, compositions of coexisting phases may be calculated and the 
miscibility gap may be plotted on a ternary diagram. However, first the compositions of the coexisting 
phases on three binary edges must be calculated. 

In a binary solution, the miscibility gap can be calculated by finding the composition of the 
coexisting phases that together represent the minimum free energy of the system. This may be done 
graphically by the tangent method, i.e., by drawing a tangent through the two points representing the 
two minima in the plot of the free energy of mixing against composition. Alternately, the relations 

^1 = ^i 

and 


may be considered. For the binary regular solution, there is a symmetric miscibility gap and therefore 


X l +X 2~ * 


X l+ X 2~^ 


x a - x@ 
X 1 X 1 


and 


2 2 


Therefore, 


ft 7 In (1 - * r ) + *2 IF = ft r In (i - * 2 ) + *2 r 

ft T In + (1 - X x ) 2 W - ft T In % 2 + (1 - * 2 ) 2 W 
Substituting x 2 = 1 - Xj into Equation 2-17, 

W 1 1 ” x i 

YT ~ 1 - 2 *, ln . 


(2-17) 

(2-18) 


(2-19) 


Equation 2-19 may be solved by an iteration method to find the miscibility gaps on the binary edges in 
a triangular diagram. 

A computer program to solve Equations 2-16 numerically and the method to form a miscibility 
gap have been presented by Kaufman and Bernstein (1970). Some examples to illustrate the possible 
solutions of certain mineralogical problems are presented elsewhere. 


14 



Chapter 3 


COMPOSITION OF COEXISTING PHASES 


IDEAL SOLUTION MODEL 


Distribution of a Component Between Two Ideal Binary Crystalline Solutions 

Although there are no strictly binary silicates, certain minerals such as orthopyroxene and olivine 
may be assumed to be quasi-binary. Because Fe 2+ and Mg 2+ are similar in ionic charge and size, 
olivine and orthopyroxene may be assumed to be binary ideal solutions. This assumption will be 
reexamined later. 

Ramberg and DeVore (1951) considered the following ion-exchange equilibrium between olivine 
and pyroxene: 

Mg Si 0 3 Si O^Fe Si C> 3 + lMg 2 Si 0 4 . (3-1) 

en fa f s f o 


The equilibrium constant for this reaction at a certain P and T is 


K 


3-1 


Fe 




( 1 - * 


op X 
Fe 


) * 


o 1 
Fe 


(3-2) 


The equilibrium constant K is a function of P and T only. In the present case, however, A 31 is not 
found to be constant except at high temperatures. (See Olsen and Bunch, 1970.) 

It may be noted that Equation 3-1 is written on a one-cation exchange basis. It may also be 
written as 


2 Mg Si 0 3 + Fe 2 Si 0 4 ^ 2 Fe Si 0 3 + Mg 2 Si 0 4 - (3-3) 

The equilibrium constant for this reaction is 


K 


3-3" 


/ opx \ 2 / 1 01 X 

U ) (1 -x ) 

F e Fe' 

1 _^opx 2 x .l 
Fe Fe 
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A Roozeboom figure with such K values has been presented by Kern and Weisbrod (1967, p. 224). It 
is known empirically from the distribution data in several mineral assemblages that equilibrium con- 
stants or distribution coefficients such as K 3 3 are very cumbersome to handle and inconsistent with 
petrological observations. One may, therefore, prefer to use the distribution data on a one-cation ex- 
change basis. It is obvious that in actual calculations of the energy values, it will be necessary to adjust 
for the activity-composition relations, such as Equations 1-7, 1-8, 1-14, and 1-16, as discussed before. 

Generally olivine and pyroxene coexist with several other minerals of fixed or variable composi- 
tion. If no significant change in the concentration of the minor components changes the binary charac- 
ter of the two minerals, ATj.j is not a function of any changes in the number or proportion or composi- 
tion of other coexisting phases. This is generally true about equilibrium constants in other systems also. 
At a certain P and T the stability of the olivine and pyroxene combination is a function of the presence 
or absence of quartz, but the value of x itself is not affected. 

Kretz (1959) used Roozeboom plots extensively to show the orderly distribution of cations be- 
tween coexisting silicate minerals in rocks. If chemical equilibrium is closely approached in the distri- 
bution of a component between two binary solutions at a certain P and T, the distribution isotherm is 
a smooth curve. If at the same time both the solutions are ideal, it will be of the form shown in 
Figure 3-1. 



Figure 3-1— Distribution of a component A between two 
ideal crystalline solutions a and 0. The numerical values are 
equilibrium constants. 
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Coexisting Ternary Ideal Solutions 

Consider two coexisting ternary phases a and P with the formulae (A, B, C)M and (A, B, C)N. The 
pure components are AM, BM, and CM in a and AN, BN, and CN in p. The chemical potentials of the 
components in a and P are 

fa=li M + RTlnx* m (3-4) 

and 

^an = P- M + R T ln *an • (3 ' 5) 


The others are similarly defined. The potentials of all the pure components /i are functions of P and T 
only, x" and x^ may be substituted forx“ M andx^ M , respectively, without altering the results. (See 
discussion before.) 


The distribution of A between a and P may be represented by the ion exchanges 

A a + B a + A jB 

and 

A a + C C a + A/3 . 

The equilibrium constants may be written as 

x a x@ 

_ B A 

* 3 - 6 ~ 


and 


A B 


X a X P 

X C *A 

K, = 3 

3-7 x a x& 
X A X C 


(3-6) 

(3-7) 


where x A = A/(A + B + C), and the other x’s are defined similarly. Both K and K 3 7 will be constants 
for all ratios of A to B to C. A plot of x A against x^ will produce a symmetric ideal distribution curve. 


NONIDEAL SOLUTIONS 


Distribution of a Component Between Two Simple Mixtures 

For the ion-exchange equation 

Aa + B/8 A (3 + B a 

at equilibrium, 

+ ^A =^B + ^A • 


(3-6) 


(3-8) 
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If (A, B)a and (A, B)j3 are simple mixtures, 


fj-l = p B “ + R T In (1 - *“) + if (%“) 2 , (3-9) 

and the other p’s are similarly defined. Substituting the values found by Equation 3-9 into Equation 
3-8 and rearranging, 


where 


4 o - *;> 


W a 
R T 


(1 



R T 



- A G° 

a 

R T 


(3-10) 


Or 


A G° - p Ba + pA 3 -p B ^ - p Aa . 


lF a IF/3 

ln ln K D - Tr o - 2 *A> + rr f 1 - 2 «5) 


(3-11) 


where = exp ( -G°/RT) and is the distribution coefficient. 

If a good least-square fit can be obtained for the distribution data by using Equation 3-1 1, it may 
be found that both minerals are close to being simple mixtures. 

As one or both of the minerals becomes less ideal, the distribution isotherms may attain different 
forms. (See Mueller, 1964.) Figure 3-2 shows an example where one of the mineral’s a is ideal and /3 
is nonideal. W 13 is assumed to vary linearly with 1/72 The values of W P /RT and K a at 673 and 1673 K 
are 2.75, 1.603, 0.77, and 1.518, respectively (see Saxena, 1969). The forms of the distribution iso- 
therms are very different from the symmetric ideal curves. 



Figure 3-2— Distribution of a component A between 
an ideal solution a and a regular solution 0 for the 
following data: 


7" (K) 

V\fi 

RT 

K 

673 

1673 

2.75 

.77 

1.60 

1.52 
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Coexisting Regular Ternary Solutions 


The composition of two coexisting phases that obey the same equation of state is considered here 
as an example. These phases are products of unmixing in a ternary solution (A, B, C)M. , W BC , 
and W' KC are assumed to be 6300, 29 000, and 38 000 J/mole (1500, 7000, and 9000 cal/mole), 
respectively , and the values of W' are assumed not to be functions of P, T, and composition (regular 
solution). Figure 3-3 shows the miscibility gap in the system and the tie lines for the coexisting phases. 
Let the phase rich in C be denoted by a. and the phase poor in C by £?. For the chemical potentials, 

+ R T \ n f a x a ( 3 - 12 ) 

and 

+ RT In fP x P .... (3-13) 


Any one of the following ion exchanges between a and ($ may be considered: 


and 


A j3 + B a^A a + B fi , 

B a + C — B /? + C a , 

Aa + C/3^Ca + A/3. 



(3-14) 


Figure 3-3— Coexisting regular ternary solutions. 
The components are A, B, and C. and 

|/V' AC are 6300, 29 000, and 38 000 J/mole (1500, 
7000, and 9000 cal/mole), respectively. The temper- 
ature is assumed to be 1573 K. 
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The equilibrium constant for Equations 3-14 is 


*A*B 


x{ x “ 

A B 


fjjj 

ft f B 


= K 


3-14 


(3-15) 


In this particular case because a and obey the same equation of state, A G° = 0 and A^ 3 . l4 = 1. In 
other cases where a and § are minerals with different crystal structures, the equilibrium constant is 
not equal to 1. The/ terms in Equation 3-15 are functions of P, T, and the ratio of A to B to C. There- 
fore K d x|) also changes with P, T, and the ratios of A to B, B to C, and A to C. 

Let the ratio of A to B to C change systematically as listed in Table 3-1 . A plot of Xg against x^ , 
where x is either A/(A + B) or A/(A + B + C), shows a smooth distribution curve (Figure 3-4). The 
form of the curve, however, is markedly different from the ideal distribution curve. 

The activity coefficients are given by 


and 


K T ln/ A 
R T lnf B 

R T In f Q 


= (* B ) 2 W 'AB + (*C > 2 ^AC + *B *C ( W '*B ~ 

- (* c ) Hgf. + (* A ) 2 ^AB + x c X A (®BC ~ ^AC 

= ( x a ) 2 ^AC + (*B ^ 2 ^BC + X Pi X B (^AC “ W *B 



+ *«)• 


(3- 16a) 
(3- 16b) 

(3- 16c) 


where x A = A/(A + B + C) and the other x’s are similarly defined. It may be checked that substitution 
of /values into Equation 3-15 gives the equilibrium constant as unity. 


Table 3-1— Composition of coexisting phases in ternary regular solutions. 


a 

P 

Y® 

X P 

k d 

B 

A 

C 

B 

A 

C 

0.010 

0.082 

0.908 

0.043 

0.868 

0.089 

0.108 

0.047 

0.406 

.040 

.079 

.881 

.174 

.718 

.109 

.336 

.195 

.478 

.080 

.071 

.849 

.339 

.522 

.140 

.530 

.394 

.576 

.100 

.066 

.834 

.413 

.431 

.156 

.602 

.489 

.632 

.130 

.056 

.814 

.512 

.309 

.179 

.699 

.624 

.715 

.170 

.037 

.793 

.625 

.171 

.205 

.821 

.785 

.796 

.200 

.020 

.780 

.700 

.081 

.219 

.909 

.896 

.862 
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Figure 3-4— Distribution of a component B between two 
ternary regular solutions plotted on a Roozeboom diagram, 
x = B/(A + B). 


DISTRIBUTION OF A CATION BETWEEN TWO OR MORE MULTICOMPONENT MINERALS 

Many rock-forming minerals are complex multicomponent crystalline solutions. The distribution 
of cations in two or more coexisting minerals in natural assemblages may still yield certain valuable 
information. The method to be followed in such cases has been discussed by Kretz (1959). In silicates 
there are at least two types of coordination for cations. Si 4+ , Al 3+ , Fe 3+ , and less commonly Ti 4+ are 
in tetrahedral coordination. Fe 2+ , Mg 2+ , Fe 3+ , Al 3+ , Mn 2+ , and Ti 4+ are found in octahedral 
coordination. Such differently coordinated ions may be regarded as forming submixtures. The distri- 
bution of Fe 2+ or Mg 2+ or any other octahedrally coordinated ion may be examined in two or more 
such submixtures forming parts of different minerals. It should be noted, however, that the chemical 
potentials of a cation in octahedral coordination may also be a function of any chemical variation in 
the concentrations of the tetrahedrally coordinated ions. Such information can be usually obtained 
beforehand by considering the chemical composition of individual minerals. For example, the positive 
correlation between the concentration of tetrahedrally coordinated Al 3+ in amphiboles and biotite 
with the Fe 2+ /Mg 2+ ratio in the mineral is now well known (Ramberg, 19526; Saxena, 1968a). 

It may be argued that the study of the distribution of a component between only two of the co- 
existing minerals that are quasi-binary solutions out of an entire assemblage of five or six minerals 
could not be useful. That is, the presence or absence of a third or fourth mineral in the assemblage 
ought to affect the distribution coefficient. This is not generally true. The distribution coefficient 
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changes only when the presence or absence of a third mineral is associated with a significant change 
in the concentration of one or more elements in one or both of the coexisting minerals. For example, 
Ti0 2 is only sparingly soluble in olivine and orthopyroxene. The chemical potential of Ti0 2 may 
increase or decrease in the rock, and rutile may be added or removed from the assemblage, but K D for 
the distribution of Fe 2+ and Mg 2+ would not change. However, if the change in M Tic , 2 changes the 
concentration of Ti0 2 significantly in one of the two coexisting minerals, K D may also change. Thus 
it is only meaningful to consider the concentrations of all the components in the two minerals and not 
the presence or absence of another phase or the change in the bulk composition of the rock. 

One of the important results of the study of cation partitioning is the recognition of how closely 
chemical equilibrium may be approached in the rocks. Whether the minerals are ideal or not, the dis- 
tribution of a component between two coexisting binary phases at a certain P and T will be repre- 
sented by a smooth distribution curve if chemical equilibrium is closely approached. If the minerals 
are not binary, the concentration of other components because of the diadochic or substitutional 
relationships may affect the orderly distribution as discussed before. In fact, the approach to chemical 
equilibrium can be studied with respect to each component individually. Figure 3-5 shows the distri- 
bution of Mn in coexisting minerals from charnockites (Saxena, 19686). Such orderly distribution of 
Mn is common in other rocks as well. The distribution of Fe 2+ and Mg 2+ between coexisting olivine 



■ bi X Capx • opx o gar 


Figure 3-5— Distribution of Mn in coexisting minerals 
in charnockites of Varberg, Sweden (Saxena, 19686). 
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and orthopyroxene at 1073 and 1 173 K was experimentally studied by Medaris (1969). Although 
Medaris made repeated grinding and heating of the reaction products, Figure 3-6 shows that the distri- 
bution points both at 1073 and 1 173 K show some scatter. The difficulties are related to the kinetics 
of the ion-exchange reaction as equilibrium is approached, particularly when the distribution approaches 
a 1 to 1 ratio in the two minerals. In contrast to these experimental results, the partitioning of Mg 2+ 
and Fe 2+ between orthopyroxene and Ca pyroxene (Capx) in metamorphic rocks as studied by Kretz 
(1963) is remarkably orderly. Most distribution points fall on a smooth curve, and a distribution curve 
representing igneous rocks is clearly separated from a distribution curve for the metamorphic rocks. 

Studies of partitioning of cations between coexisting minerals in natural rocks by petrologists 
(Albee, 1965; Annersten, 1968;Binns, 1962; Butler, 1969; Kretz, 1959; Gorbatschev, 1 969; Hietnan, 
1971 ; and Mueller, 1960, among many others) are attempts to rationalize the concept of metamorphic 
facies on a mineral and chemical basis. In experimental systems, similar attempts have been made by 
Nafziger and Muan (1967), Larimer (1968), and Medaris (1969), among others. The results of such 
partitioning studies have generally confirmed the usefulness of the approach and the need for more 
thermodynamic data on crystalline solutions. 

In essence, problems of phase equilibria are distribution problems, and a statistical approach to 
such problems may be made to avoid the consideration of the thermodynamic properties of solutions 
in individual minerals. Such approaches have been made principally by Greenwood (1967) and Perry 1 
and should be applicable in solving the petrogenetic problem of incompatible assemblages and the 
recognition of chemical equilibrium in natural or experimental systems. 
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Figure 3-6— Distribution of Fe 2+ and Mg 2+ 
between synthetic olivine and pyroxene. Distri- 
bution data from Medaris (1969). 


1 


K. Perry, Jr.: “Construction of a Single (m + 2) Dimensional Phase 


Diagram From Petrochemical Data.” 


To be published. 
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Chapter 4 


MEASUREMENT OF COMPONENT ACTIVITIES BY ANALYSIS OF PHASE DIAGRAMS 


An experimental measurement of activities of components in a crystalline solution, particularly the 
silicates, is beset with difficulties and the measured values are subject to large errors. Therefore, obtain- 
ing such activity-composition relations from phase diagrams would be very convenient. 

However, there is no direct method of doing this without some kind of a solution model. The use 
of a solution model brings in uncertainties in the activity values, which depend in extent and form on 
the choice of the model. The attempts to obtain the thermodynamic functions of mixing through the 
use of various solution models is still useful. For some crystalline solutions it may be possible to 
experimentally determine such properties. A comparison between the experimental values and the 
values based on a particular model would provide a greater understanding of the interrelationship of 
the crystal structural parameters on which the model is based and the thermodynamics of the crystal 
phase. In other cases where experimental determinations cannot be made, the empirically derived 
functions of mixing may be tested for their physical significance by their success in application to prob- 
lems of petrogenesis. 

SYMMETRICAL MIXTURES 

Simple Mixture 

The thermodynamics of a simple mixture or regular symmetric solution has been discussed before. 
Assume that the solution (A, B)M unmixes into two coexisting solutions, a, rich in AM, and |3, rich in 
BM. At equilibrium, 

and (4-1) 

so that according to the simple mixture model, 

+ R T In *“ + W (1 - *“) 2 =/iJ> + R T In + W (1 - x^) 2 (4-2) 

where and (3 Q stand for the same pure end member structure AM. Eliminating and and sub- 
stituting jc P A = 1 - x ^ , 

R T in ** + W (1 - *“) 2 = R T In (1 - *“) + W (*“) 2 (4-3) 
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or 


(4-4) 


W 

RT 


_ ^ [(1 -*“)/*“] 


2 *A 


This expression is similar to the one obtained by Thompson (1967). The equation for the curve of co- 
existence of two phases may also be written in terms of critical temperature T c of unmixing and the 
mole fractions by substituting 

W = 2 R T 

C 


into Equation 4-4: 


T = 2 T 


1 - 2 ** 


In [(!-*£)/*;] 


(4-5) 


If there are data on the composition of coexisting phases at different temperatures and the form 
of the solvus is symmetric, the value of W and the activity-composition relations can be calculated. 


Symmetrical Mixture of Higher Order 

Symmetrical crystalline solutions may not be simple mixtures and may require an expression with 
two or more constants to represent G £M : 

= o-*;> [a, + (i - 2 . <«> 

For such a solution, Equation 4-3 is 

R T In ** + [Ao + ^ (1 - 2 **) 2 ] (1 - **) 2 = R T 1 n + [d Q + T, ( 1 - 2 *£) 2 ] (1 - * J)2. 

(4-7) 

Using the relation p B = ju^ , an equation similar to Equation 4-7 can be written, and the two equations 
can be solved simultaneously to obtain A Q , A 2 , and the activity-composition relation. 

ASYMMETRICAL SOLUTIONS 


J 


Subregular Model 

As mentioned before, may be expressed as a polynomial in the mole fraction x A or x B for 
the compound (A, B)M according to Guggenheim’s equation: 

G BM = *A ( X A - * B ) + ^2 - V 2 + • • •] • (1 ' 49) 

If A 2 = 0 is substituted into Equation 1-49, a two-constant equation for an asymmetrical solution is 
the result. Proceeding as in the previous sections, 

R T In *“ + R T In /"* = R T In x? + R T In f? (4-8) 
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and 


R T In x* + R T In f* = R T In x? + R T In fP. (4-9) 

Substituting values of R T In / from Equations 1 -5 0 and 1-51, 

KT In + (*£> 2 [4 0 + 4, (3 ** - *£>] = K T In + (,^) 2 p 0 + A t (3*J-*g)] (4-10) 

and 

R T ln *b + ("a) 2 Po - ^ ( 3 *b ~ *a)] = « T ^ 4 + (^) 2 K - ^ (3 - *J)] • (4*1 1) 

Equations 4-10 and 4-1 1 now can be solved simultaneously to yield the values of the two con- 
stants A Q andA l . 

The method of calculation presented above is equivalent to that used by Thompson (1967) and 
Thompson and Waldbaum (1969a, b). Thompson’s (1967) equation for G £M is 


^EM ~ *A ^2 + *B ^1 


(4-12) 


where 


and 


(j i x a x b w g~ 


( 'l x A X B • 


This is as if the crystalline solution is composed of x A moles of a simple mixture with IVg 2 and x } 
moles of another simple mixture with W G . Then, 

C EM = ( W G X X B + W G 2 *A> X A X B 


[ v i -‘*> + v ! -* b >] 


X. X TA 

A B 


2 - 2 X, 


+ K, 


2-2 


1 2 2 


X X„. 
A B 


(4-13) 


Substituting 1 =x A +x B into Equation 4-13, 


C EM ~ 


' v G W G 

~Y C 1 - *A + *B> + — 1 t 1 + *A - V 


X. x n 

A B 


W G + w w _ w 

°2 + G 1 G 2 G 1 . 

+ (x. - x„) 

n T o 'A B' 


X . *_ 

A B 


(4-14) 


which is of the same form as Guggenheim’s equation with two constants A Q and A j . 
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Therefore, 


and 


=■ 


w o t + 


(4-1 5a) 



(4-1 5b) 


A 0 /RT and A l /RT would correspond to the notation B G and C G , respectively, used by Thompson 
(1967) following Scatchard and Hamer (1935). 

For the activity coefficient, 

RT \n f k - ( x B ) 2 U Q + A l (3 * A - * B ) + • • *3 • (1-50) 

Substitution of Equations 4-15 into Equation 1-50 gives 


*Tlnf A = (x b )2 


"g 2+ % 


( x A ~ X B + 2 **) 


= x 2 

B 


w r + w r w r _ r, 

G 2 t Gj G 2 Gj 


<*A - V 


+ ( W G 2 “ ^Gj) * A | 


= *1 t(» e , *B + % *A> + (*<,, ‘a-'g,**)] 

= *B »G, <*B - «A> + 2 »G, *A> 

- *b , n - 2 * A ) + 2 ip G2 * A j 

= *B »G, + 2 -A - * 4 " 16 ) 


which is the same form used by Thompson (1967). 

Substituting the calculated values of A Q and A j into Equation 1-49, ^EM can be estimated. The 
solvus bounding the two-phase region can then be determined by the graphical tangent method or by a 
suitable iteration numerical method. The calculated values of x B> , and x| are then compared 
to the observed mole fractions to test the applicability of the model. Other excess functions of mix- 
ing may be calculated by Equations 1-52 and 1-53. 


Quasi-Chemical Approximation 

This model has been used by Green (1970) to study the halite-sylvite solvus. Consider two coexist- 
ing phases M and N with components A and B. At equilibrium at a certain P and T, 
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and 


Ma = A 


or 


and 


U M = „N 


+ R T In *“ + R T In = /a AN + /? T ln *n + R T ln yN 
+ R T ln *“+fin n /“,uBN +ftnni N +sr ln yN . 

d B ' B B 


(4-17) 


Because both M and N obey the same equation of state, the chemical potentials of pure A in M and A 
in N and B in M and B in N are canceled. 

Substituting values of and from Equations 1-35 and 1-36, respectively, into Equations 4-17* 


2 

ln**[ + - ln 


4>l (fi- 1) 

z 

— In r N + In 

~ ^ (F - 1) 

0 “ 0 s + !)_ 

■ A 2 

(F + i)_ 


and 


(4-18) 


2 

In + — — 

B 2 


In 

r <p k (P - i)i 

1 + A 

z 9 1 

= ln + ln 

r <t>i - i)i 

i + a 


<t>\ (J3 + 1)_ 

B T — 

+ !)_ 


(4-19) 


where 0 artd /?', which correspond to phases M and N, respectively, and the 0’s are defined by Equations 
1-30 (usihg 0 A instead Of x A , etcO and 1-37. q l and q 2 are the contact factors discussed before. They 
are not independent and should approach 1 simultaneously. Green (1970) assumed \/ q x q 2 = 1. The 
t\ta> indepehdent reiatidns 4-18 and 4-19 contain two unknowns q 1 /q 2 and W and can be solved by an 
iteration process. The ratio q ] lq 2 is a function of the geometry of the substituting chemical species 
and therefore ihay be regarded as almost independent of T. Substitution of q l /q 2 back into Equations 
4-18 and 4 1 1 9 gives two independent values of W at each temperature. Any difference noted in the 
tWO Values of W would be caused by the inadequacy of the solution model to fit to the experimental 
data* 

The sOlVus bounding the two-phase region may be determined graphically by the double tangent 
ihethod on a plot of free energy Of mixing against rnole fraction or by a numerical iteration method. 
The dxcess functions of mixing can be calculated by Equations 1-38 and 1-39. 


EXAMPLE OF CALCULATION OF FUNCTIONS OF MIXING: THE C a W0 4 -SrW0 4 SYSTEM 

Chang (1967) presented the data on the two-phase regions with solvus in the binary tungstate 
R n W0 4 type of crystalline solution. Table 4-1 shows the data on the composition of the coexisting 
phases a and 0 rich in CaW0 4 and SrW0 4 , respectively. Calculated values of A Q /RT and A^/RT 
according to the stibregular model are listed in Table 4-2. The relationship between A Q /A j and T is 
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Table 4-1— Chemical composition of unmixed phases in the 
system CaW0 4 -SrW0 4 , a 


T (K) 

X Si 

y.0C 

x Ci 

l-H 

X 

x 0 

X Ca. 

823 

0.005 

0.995 

0.995 

0.005 

873 

.010 

.990 

.980 

;020 

923 

.025 

.975 

.955 

.045 

973 

.035 

.965 

.905 

.095 

1023 

.067 

.933 

.800 

.200 

1073 

.120 

.880 

.630 

.370 


a a and /3 are coexisting phases rich in CaW0 4 and SrW0 4 , respectively. The 
compositions are from Chang (1967, Figure 3). 


Table 4-2— The calculated A 0 /RT and 
AJRT in (Ca, Sr)W0 4 . 


T (K) 

Aq/RT 


600 

4.364 

-0.327 

650 

3.524 

-.348 

700 

3.119 

-.455 

750 

2.537 

-.514 

800 

2.021 

-.607 

900a 

0.839 

-.739 

1000 a 

-0.296 

-.884 


a From Equations 4-20 and 4-21. 

Note: An error of ±5 percent in the mole fractions 
(Table 4-1) results in a ±800 J/mole (±200 cal/mole) 
error in determining A 0 and A j . 


linear and is given by 

14.1526 - 0.011 35 T 

0.9616 - 0.001 45 T . 
can then be calculated from the relation 

^EM ~ *SrW0 4 *CaW0 4 ^0 + ^1 (*SrWC> 4 “ *CaW0 4 ^’ 

and the activity coefficients from the relations 

^ ^ f SrW0 4 = *CaW0 4 ^0 + *SrW0 4 ~ *CaW0 4 ^ 


and 


l o 


RT 
A 


l 


RT 


(4-20) 

(4-21) 


30 






and 


^ ^ In ^CaW0 4 *SrW0 4 ^0 ^1 ^ *CaW0 4 “ ^SrWO^ ' 

Figure 4-1 shows the activity-composition relation at 1073 and 1273 K. G EM can also be plotted 
against composition, and the composition of the coexisting phases can be found by the tangent 
method. In the present case, the differences between compositions calculated by the model and the 
observed compositions in Table 4-1 are found to be small. 

As mentioned before, the activity-composition data and other thermodynamic functions as cal- 
culated from phase diagrams are sensitive to the nature of the assumptions and the model used. For 
the system CaW0 4 -SrW0 4 , calculation shows that the use of quasi-chemical approximation predicts 
the solvus with the same accuracy that the subregular model does. The use of quasi-chemical approxi- 
mation requires the values of , q 2 , and the coordination number z, the number of the nearest Ca 2+ 
or Sr 2+ ions surrounding each other. In CaW0 4 there are four Ca 2+ ions surrounding each Ca 2+ at a 
distance of approximately 0.39 nm. There are four more Ca 2+ ions at a distance of 0.5 nm. z may be 
assumed equal to 4, and the two equations of the quasi-chemical approximation (Equations 4-18 and 
4-19) can be solved simultaneously to find qjq 2 . 

In this case, it may be assumed that either q x + q 2 = 2 or s /q l q 2 = 1. The differences in the cal- 
culations of W using either of the two assumptions are small. (See Green, 1970.) A computer pro- 
gram may be used to solve each of the two equations independently by using various values for q x and 
q 2 and to compare the W values so obtained until the best set of W values is found. The W values with 
z = 4 are listed in Table 4-3. q x and q 2 are 1.20 and 0.80, corresponding to SrW0 4 and CaW0 4 , 
respectively. The atomic radii (R) for Sr 2+ and Ca 2+ are 1.12 and 0.99, respectively (Ahrens, 1952); 
therefore, the ratio of q l to q 2 is 1.50 and the ratio of the radii of Sr 2+ to Ca 2+ is 1.13; these values 
are not similar. 



Figure 4-1— Activity-composition relation in CaW0 4 - 
SrW0 4 crystalline solution at 1073 and 1273 K. 
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Table 4-3—2 W/zRT for the system 
CaW0 4 -SrW0 4 . 


T (K) 

Equation 

4-18 

Equation 

4-19 

873 

5.883 

mm 

923 

5.173 

Hi 

973 

4.919 

4.460 

1023 

4.495 

4.301 

1073 

4.183 

4.386 


z = 4;q j = 1.20 ;? 2 = 0.80. 


The following equation describes the relation between the calculated W (the average of the two 
values listed in Table 4-3) and T: 

^ = 34.64 - 0.055 06 T + 0.000 025 T 2 . ( 4 - 22 ) 

Figure 4-2 shows a comparison of G EM at 1073 K calculated according to both the subregular 
(SR) model and the quasi-chemical (QC) model. The value of G EM according to the latter is nearly 
twice that calculated according to the former. Differences between the other calculated functions of 
mixing, H EU and *EM , are even more marked. Unfortunately there are no data on experimentally 
determined // EM and for the CaW0 4 -SrW0 4 system, and, therefore, there is no way to know 
which model predicts the thermodynamic functions better in this particular case. 

For the system NaCl-KCl, Green (1970) compared the thermodynamic quantities calculated by 
the subregular model and the quasi-chemical model with those determined by experiments. The 
thermodynamic quantities predicted by the quasi-chemical model are closer to those measured 
experimentally. 

A comparison of the predictions of the functions of mixing in several binary alloys by regular 
solution model and by the quasi-chemical model (Lupis and Elliott, 1967) shows that generally the pre- 
dictions by the latter for the excess free energy are closer in agreement with experimental determina- 
tions than those by the former. The prediction for the excess entropy from the quasi-chemical model 
is not satisfactory. This may be in part caused by the neglect of the nonconfigurational excess entropy 
in many of the binary alloys. For the halite-sylvite system, Green (1970) finds that the nonconfigura- 
tional contributions are unimportant and suggests that the positive excess entropy of mixing found in 
the NaCl-KCl system may result from the introduction of vacancies or other defects into a crystalline 
solution. 

This approach of calculating thermodynamic functions of mixing by the analysis of phase diagrams 
is relatively new in the field of mineralogy and deserves more attention from mineralogists and petrolo- 
gists. The fact that there is no unique analysis of a solvus and two or more solution models may be 
applicable to the same solvus data need not be a barrier to acquiring and interpreting more phase data 
with the help of various solution models. Experimental verification of many of these results may not 
be possible in the near future. However, it may be possible to test such thermodynamic data by their 
application to petrogenetic problems and get physically meaningful results. 
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Figure ^-Thermodynamic excess functions 
of mixing in (Ca, Sr)W0 4 according to the 
two-constant asymmetric model and the 
quasi-chemical approximation. 
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Chapter 5 

MEASUREMENT OF COMPONENT ACTIVITY USING COMPOSITION 
OF COEXISTING MINERALS 


Experimental data on the distribution of a component between two coexisting crystalline solutions 
at a fixed P and T for systems such as olivine and pyroxene have been collected by Nafziger and Muan 
(1967), Larimer (1968), and Medaris (1969). Distribution data are also available for natural assemblages, 
but the P and T of their formation are indefinite. The distribution data from natural assemblages in 
many cases may be found to represent ion-exchange equilibrium closely. If precise P and T are not 
important, such data may be used to obtain useful information on the thermodynamic nature of mixing 
in the minerals. For this purpose, the thermodynamic equations according to various solution models 
for binary solutions presented in this section may be used. 

The composition of coexisting phases that do not obey the same equations of state may be used 
to find the activity-composition relations in each phase in suitable cases. Consider a and j3 with chemi- 
cal formulas (A, B)M and (A, B)N, respectively, which are in ion-exchange equilibrium at a certain P 
and T. From Chapter 3, 

A a + B B a + A/3.. (3-6) 


The equilibrium constant K 3 ^ is given by 

„ _(< <\fii 
3 \-i4Nifi) 


(5-1) 


The term in the first bracket is the distribution coefficient K D . Depending on the nature of the data 
available, the following cases may be considered. 


COMPOSITIONAL DATA AVAILABLE ON A COMPLETE DISTRIBUTION ISOTHERM 


The simple mixture model, the two-constant asymmetric model, and the regular solution model 
with quasi-chemical approximation are the possible choices. According to simple mixture model, 

W a n 

In K = In A - _(1 - 2 jO + — (1 “ 2 *£). (3-11) 

3-6 D ft y A' R T ^ A 

A nonlinear least-square fit using the data on x^ and x^ finally yields K 3 _g , W a , and W 13 . 
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According to the Redlich and Kister equations (King, 1 969, p. 326), 



R T ln f A = X B l A 0 + A 1 ( 3 X A ~ *B> + A 2 ( X A ~ X B> ( 5 X A ~ *B> + ' 


(5-2) 

and 





R T In f h = x\ [a o - A 1 (3 * R - * A ) + A 2 (* b - x A ) (5 * B - * A ) + • 


(5-3) 

Therefore, 





, / A 

R T ln__= A 0 (* B - * A ) + (6 x A * B - 1) + A 2 (% B - * A ) (1 - 8 * A * B )- (5-4) 

' B 


Substituting the values of f\tf\ and f 0 /fg by using Equation 5-4 in Equation 5-1 , neglecting the 
constants A , ’s, arid rearranging in logarithmic form yields 


A a 

In K = In K + — — (*“ - 
3-6 D R T A 


O + 


A\ 

— (6 

R T K 


x a x a 
X B *A 


A? 

n o 






/? T 


R T 


(5-5) 

Equation 5-5 is of the form 

IV = M + 4 + A 2 * 2 + /I 3 x 3 + . . . 

where the TV and the x’s are known quantities, It may be solved by a numeric least-square method 
yielding M- In and other constants. There must be a minimum of five distribution points. 

According to the quasi-chemical approximation, 


f.= 


1 + 


0 B (P - .!) 

4 >„ 0 + 1)1 


z< ? a /2 


(1-35) 


and/ B is defined as in Equation 1-36. Substituting these values of/ into Equation 5-1, 


\ ^ <t>l 0 a - 1)' 

Z X /2 


* P qy2 

4>l 0 a + 1) 


• 0 P + l) 


4>l 0 a - i) 

1 + 

z a q A / 2 

", t K 0 8 * - !)' 


(/? a + i) 


+ 0g (/3* + 1). 



where (5° and $ are for phases a and 0, respectively, and are given by Equation 1-30, with x A replaced 
by <A a , etc. q a and q 0 are contact factors for phases a and 0, respectively. A numerical least-square 
method may be used to solve Equation 5-6. 
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COMPOSITION DATA ON A COMPLETE DISTRIBUTION ISOTHERM AND THE 
ACTIVITY-COMPOSITION RELATION IN ONE OF THE TWO COEXISTING PHASES 


Depending on the accuracy of the data, Equations 3-1 1, 5-5, and 5-6 may be used. If necessary, 
all three constants in Equation 5-4 may be used. Equation 5-4 may be written as 

(2A 0 +6A 1 + lO A 2 ) - x\ (6 A x + 24 4 2 )+% 3 ( 16 A 2 ) - (A Q + A x + A 2 ) . (5-7) 

’ B 


Transforming Equation 5-1 into logarithmic form and substituting values from Equation 5-7 for 


fiifl • 


t a 


A n + A. + A , 


lnK n+ ln_ = lnL, + J! I 


f a 

1 A 


3-6 


R T 


2 A 0 + 6 A x + 10 AA 

rJ J 


16 A 2 
~R T 


(5-8) 


Equation 5-8 is of the form 


Y = a Q + AjX + a 2 x 2 + a 3 x 3 + 


If the activity-composition relation in a is known, Equation 5-8 may be solved by least-square analysis. 

If sufficient number of distribution points are available, Equation 5-5 may be used with the third 
constant A 2 and the results compared with those obtained by using Equation 5-8. 

The olivine-chloride solution system (Schulien et al., 1 970) has been used by Saxena 1 for calcu- 
lating the activity-composition relation in binary solutions using the distribution data. 


\ yllfir.oHiv'i'’- 



*S. K. Saxena: “Retrieval of Thermodynamic Data From a Study of Inter-Crystalline and Intra-Crystalline Ion Exchange 
Equilibria.” To be published in Amer. Mineral., 1972. 
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Chapter 6 


ORDER-DISORDER IN Fe 2+ -Mg 2+ SILICATES 


Long-range order-disorder phenomena in ferromagnesian silicates differ from those in alloys in 
several important respects. First, as opposed to alloys, the silicate framework remains more or less 
inert and only a certain number of cations take part in the site exchange. Second, a complete crystal- 
line solution usually exists between the Mg and Fe end members. Third, because Fe 2+ and Mg 2+ are 
similar in size, charge, and other characteristics, the site preference energies (corresponding to the dif- 
ference in binding energy of the ion between the nonequivalent sites) are not strongly dependent on 
the degree of order as is usual in many binary alloys. Order-disorder or the intracrystalline cation dis- 
tribution in silicates is measurable by X-ray (see Ghose, 1961) and other spectroscopic techniques. The 
energy of the intracrystalline ion exchange is part of the Gibbs free energy of the crystal and is, there- 
fore, a very useful thermodynamic quantity (see Mueller, 1 969, and Thompson, 1 969). 


INTRACRYSTALLINE ION EXCHANGE AND SITE ACTIVITIES 

A crystalline solution (A^ , B l x )M may have the two cations A and B distributed between two 
nonequivalent sites a and (3. M is the inert silicate framework. Following Dienes (1955) and Mueller 
(1960, 1962, 1969), the disordering process may be represented by the following exchange reaction: 

Aa + B0 ^ AJ3 + Ba . (3-6) 


In terms of kinetic theory, the time rate of change of A in site j3 is given by 


dX A 
d t 


- K /3a ^/3a '‘A ■* B 


*1 < ~ <*>,« *5 


V ^a/3 *B *A ’ 


( 6 - 1 ) 


where x refers to the mole fractions, K 0a and K a0 are rate constants and functions of P and T only, 
and are <t> 0a and <j> a0 analogous to activity coefficient products in a macroscopic chemical system and 
are functions of P, T, and 0. At equilibrium, 


d xQ 

A 

d t 


= 0 


and 


K 


_ 4 n *b n 


a 1 ® o“ 

A B 


3-6 


K 


( 6 - 2 ) 


a/3 




39 



where / is the partial activity coefficient and a the partial activity. The product of the /’ s appears as <p 
in Equation 6-1. The term “partial” is used to distinguish between the activity of A on the site from the 
activity of A in the crystal. 


The distribution coefficient is 


K. 


A B 




X 


a 

A 


K d has sometimes been referred to as the ordering parameter. The distribution coefficient, however, 
should not be confused with the ordering parameter S used to describe ordering in alloys. S' = 1 cor- 
responds to the highest possible order, and S = 0, to complete disorder. This is the opposite in the case 
of K d . Further, K D will be used to describe order-disorder in nonstoichiometric silicates forming com- 
plete crystalline solution series. In such silicates the formation of a fully ordered or disordered periodic 
structure is not possible. Even with the greatest tendency towards ordering, some of the excess atoms 
of one component must inevitably occupy sites belonging to the other, which leads to a lower order 
or disorder. The distribution coefficient is a function of T and the varying ratio of A to B in the crystal 
and, therefore, is of little thermodynamic significance. 

The equilibrium constant K 3 ^ is a function of P and T only. However, as the volume changes in- 
volved in the ion exchange are negligible, the dependence of K 3 ^ on P is ignored, and K 3 ^ is considered 
to be only temperature dependent. 

The definition of the chemical potential of a cation on a site presents certain problems (Mueller, 
Ghose, and Saxena, 1970). One may write the following equations, as done by Grover and Orville 
(1969), for the chemical potential of a cation A on the sites a and /3: 

^ = (6-3) 

and 

= ^A 0+ R Tln a B’ (6-4) 

where the m° ’s are the standard chemical potentials and the a’s are the corresponding partial activities. 
According to classical thermodynamics, however, it is incongruent to define two different chemical 
potentials for one species in a single homogeneous phase. In such a case, 

and (6-5) 

^A° 


To avoid this difficulty, Borghese (1967) regards A in site a as a distinct species from A in site /3. 
This is somewhat analogous to speaking of the chemical potentials of 0 2 and 0 3 in a homogeneous 
gas phase. The idea of defining a new potential analogous to chemical potential called a site preference 
potential (Greenwood in Grover and Orville, 1969) could also be considered. 
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The use of a site preference potential may be avoided in practice, particularly because its quan- 
titative use in thermodynamics is rarely of importance. 

The standard site preference energy or the intracrystalline ion-exchange energy A G° for the ex- 
change reaction (Equation 3-6) is given by 

A G° = - R T In K a (6-6) 

where K g is the equilibrium constant and is a function of T only, unlike K D , which is a function of 
both T and composition. 

THERMODYNAMIC FUNCTIONS OF MIXING 

One of the principal aims of the study of order-disorder phenomena is to investigate the thermo- 
dynamic properties of the crystalline solution as a whole. In case of an ideal macrophase, the activity- 
composition relation is given by 

«!=(•*-,)", (6-7) 

where N is the number of structural sites in the crystal. When there are two sites, 

a i = ( *i ) > 

and if these sites are different, 


«i = (*? + *f) 

(6-8) 

- (*i) a 

(6-9) 


where a and ft are two nonequivalent structural sites. The latter method has been generally used 
(Mueller, 1962; Thompson, 1969). Extending the above method to the nonideal case, 

a. = (a.) a (a.y, (6-10) 

where a“ and a? are partial activities referring to the sites. In an orthopyroxene (MgMgSi 2 O fi - 
FeFeSi 2 0 6 ) where there are two sites, Ml and M2, the activity of Fe 2+ in the crystal may be expressed 
as 


a opx = a M1 a M2 
2 + 2 + 2 + 
Fe Fe Fe^ 


If the activity is considered on a one-cation basis, i.e., for the crystal (MgSi0 3 -FeSi0 3 ), 

(a* 


0 °px 
Fe 2 + 


,M1 „M2 \l/2 . 

^ 2 + _ 2+2 
Fe Fe 


( 6 - 11 ) 


( 6 - 12 ) 


The partial activity af is equal toffx? where / is the partial activity coefficient. At a certain tempera- 
ture the atomic ratio x i in the two sites a and (3 can be determined by X-ray or other resonance tech- 
niques. The next problem involves the evaluation of the partial activity coefficients. 
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Several crystals of suitable composition (A, B)M between the end members AM and BM may be 
chosen and heated at a certain temperature for a time long enough to attain equilibrium for the intra- 
crystalline ion exchange (Equation 3-6). Several such distribution isotherms may be obtained. A 
model suitable for interrelating the partial activity coefficient with the atomic fraction at the site 
must be found. The simple mixture or regular solution model may be found useful in cases where the 
form of the distribution isotherms does not indicate too much of a nonideal state of mixing A and B 
at a and 0. Thus, 

In K D - 1^(1- 2**) t f^(l - 2 ,*> . (6-13) 

where W is related to the partial activity coefficient by 

R T In = W (1 - *“) 2 . (1-27) 


At this point, certain other partial functions of mixing may be considered. The partial free energy 
of mixing at the sites is given by 

G^^nin^+^nina* (6-14) 

= *A R T In /a + R T In /£ + T (*a R ln *A + x b R ln *b> * (6-15) 

Substituting ,9“ = -/?(*“ ln.v“ + x“ lnx“) and Equation 1-27 into Equation 6-15, 

G£=l (6-16) 

The term W a x °^ is also the partial excess free energy of mixing. 

In a crystalline solution such as AAM-BBM, the total free energy of mixing is given by 

X 1 + x'f *2 + 

G ~ — R T ln oj-af + — - R T ln a*- a? , (6-17) 

M2 aa t 2 bb 

which can be shown to be 


x a - 

G * = - A R Tin * 3 . 6 - r(5“ +Sfi) + W a + WfixPxP. 


(6-18) 

(6-19) 


This expression is similar to the one derived by Grover and Orville (1969) for ideal mixing at the sites. 
Note that in the above expression K 3 _g is the equilibrium constant and not the distribution coefficient, 
as in the case of ideal mixing. 

Substituting AG° =-RT ln K 3 _g into Equation 6-19, 


:ac" - r(ss + sft + (c« + rf„ ), 


( 6 - 20 ) 


42 



Thus the free energy of mixing in the crystal as a whole is a result of the energy due to the distri- 
bution of the cation between a and j3 sites, the entropy change due to the distribution of A and B 
within a and j3 sites, and finally the excess energy of mixing that is the result of the nonideal state of 
the solution in a and j3. 


The partial excess free energies of mixing at the sites are 


and 


n a _ u a c a 

°EM “ "EM °EH 

= — T S gM 


( 6 - 21 ) 

( 6 - 22 ) 


Substituting Equations 6-21 and 6-22 into Equation 6-20, 

2 ^ G + + ^ EM ” ^ (^IM + ^EM + ^IM + ^ EM ) ‘ 


(6-23) 


The thermodynamic relations presented in this and earlier sections have been used in analyz- 
ing the data on site occupancies in orthopyroxene (Saxena; 1 see also Saxena and Ghose, 1971). 


KINETICS OF ORDER-DISORDER 

Virgo and Hafner (1969) made the important observation that there is an apparent cutoff or tran- 
sition region on the temperature scale below which no more ordering or disordering occurs. This tran- 
sition temperature in orthopyroxene was estimated to be approximately 750 K. Above this temperature 
the activation energy required for diffusion to start in the direction of disordering is of the order of 
80 kJ (20 kilocalories) (Virgo and Hafner, 1969). Below this temperature the activation energy should 
be very high. This is confirmed from the measurement of order-disorder in metamorphic pyroxenes 
that cooled slowly through geological time. Figure 6-1 shows the data on the K D values for the distri- 
bution of Fe 2+ and Mg 2+ between Ml and M2 sites in metamorphic orthopyroxene. From these data 
it may be noted that no orthopyroxene shows a degree of order representing temperatures lower than 
723 K. 

Mueller (1970) proposed a two-step mechanism for order-disorder kinetics in silicates. This in- 
volves a low-temperature process with high activation energy and a high-temperature process with a 
lower activation energy. This mechanism may be responsible for ordering characteristics distinguishing 
metamorphic, igneous plutonic, and volcanic pyroxenes. The intracrystalline ion-exchange equilibria 
in igneous plutonic rocks is not ordinarily quenched at any temperature because of slow cooling rate. 
The same applies to such equilibria in metamorphic rocks. However, in metamorphic rocks attainment 
of such equilibria is possible below the transition temperature if crystallization or recrystallization 
occurs at these temperatures. Because of rapid cooling in volcanic rocks, it is very possible that the 
intracrystalline equilibria are quenched and the temperature indicated by order-disorder is not very 
much lower than the original temperature of crystallization. 


l S. K. Saxena: “Retrieval of Thermodynamic Data From a Study of Inter-Crystalline and Intra-Crystalline Ion Exchange 
Equilibria.” To be published in Amer. Mineral., 1972. 
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Figure 6-1— Temperature estimate of the ion-exchange equi- 
librium between sites in natural orthopyroxenes; 723 K 
appears to be the transition temperature below which no 
more ordering takes place because of a potential barrier. 
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